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Abstract 



The HEIDELBERG-MOSCOW experiment, which is the most sensitive double 
beta decay experiment since ten years has been regularly continued until end of 
November 2003. An analysis of the data has been performed already until May 
. 20, 2003. The experiment yields now, on a 4c7 level, evidence for lepton number 

violation and proves that the neutrino is a Majorana particle. It further shows that 
neutrino masses are degenerate. In addition it puts several stringent constraints 
on other physics beyond the Standard Model. Among others it opens the door 
to test various super symmetric theory scenarios, for example it gives the sharpest 
limit on the parameter X'm in the R-parity violating part of the superpotential, and 
gives information on the splitting of the sneutrino-antisneutrino system. The result 
from the HEIDELBERG-MOSCOW experiment is consistent with recent results 
from CMB investigations, with high energy cosmic rays, with the result from the 
g-2 experiment and with recent theoretical work. It is indirectly supported by the 
analysis of other Ge double beta experiments. Recent criticism of various kind has 
been shown to be wrong, among others by measurements performed in 2003 with 
a "^^^Bi source (^^^i?a), by simulation of the background in the range of Q/3/3 by 
GEANT4, and by deeper investigation of statistical features such as sensitivity of 
peak search, and relevance of width of window of analysis. 



1 Introduction 

Double beta decay is the most sensitive probe to test lepton number conservation. Further 
it seems to be the only way to decide about the Dirac or Majorana nature of the neutrino. 
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Double beta decay can contribute decisively to the field of neutrino physics also by 
setting an absolute scale to neutrino masses, which cannot be observed from neutrino 
oscillation experiments. 

The observable of double beta decay is the effective neutrino mass 
(m) = I E^ei^il = \^eJ\ + e^'^2|m(2)| + e'^^\m^J^J\, 
with t/g^ denoting elements of the neutrino mixing matrix, neutrino mass eigenstates, 
and (pi relative Majorana CP phases. It can be written in terms of oscillation parameters 



'1^ 







\mfj\ 


= \uj 




= \uj 



Amli + ml, (2) 



'Amia + Amii + mf. (3) 

. the half-life for Oiy/3/3 d( 
(m^)^, and for the limit on deducible in an experiment we have 



The effective mass (m) is related with the half-life for 0z//3/? decay via [Ty2 
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Here a is the isotopical abundance of the emitter; M is the active detector mass; t 
is the measuring time; AE is the energy resolution; B is the background count rate and 
e is the efficiency for detecting a signal. Determination of the effective mass fixes the 
absolute scale of the neutrino mass spectrum [HI HHj . 

The HEIDELBERG-MOSCOW experiment has been regularly continued in 2003. It 
had to be stopped, on November 30, 2003, according to contract. Unfortunately the 
Kurchatov institute did not agree to prolong the contract. The experiment is already 
since 2001 operated only by the Heidelberg group, which also performed the analysis of 
the experiment from its very beginning. 

The experiment is since ten years now the most sensitive double beta experiment 
worldwide. In this report we will describe in section II the evidence for neutrinoless double 
beta decay (Oz//3/5), found by an analysis of the HEIDELBERG-MOSCOW experiment 
including the three more years of data taking. 

The result derived from the full data taken until May 20, 2003 is 

Tl'j^ = (0.69 - 4.18) X lO^V (99.73%c./.) (5) 

with best value of — ^-^^ ^ ^0^^ V- Thus double beta decay is the slowest nuclear 
decay process observed until now in nature. Assuming the neutrino mass mechanism to 
dominate the decay amplitude, we deduce 

(m^) = (0.24 - 0.58) el/ (99.73%c./.), (6) 

with best value of 0.44 eV. This value we obtained using the nuclear matrix element of 
|25j . Allowing for an uncertainty of ±50% of the matrix elements (see |S1E]), this range 
widens to 

(m^) = (0.1 -0.9) el/ (7) 
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The result (2) and (3) determines the neutrino mass scenario to be degenerate jT^ I^. 
The common mass eigenvalue follows then to be rricom = (0.14 — 3.6) eV (99.73%). 

The new results with three more years of statistics confirm our earlier results [U El 
m El El on a higher confidence level. The signal is now seen on a 4.2a level (see section 
2). 

If we allow for other mechanisms (see |T7l|THllinilini), the value given in eq. ©,(0) has 

to be considered as an upper limit. In that case very stringent limits arise for some other 
fields of beyond standard model physics. To give an example, it gives the sharpest limit 
on the Yukawa coupling X[^i in the R-parity violating part of the superpotential [2^. It 
also gives information on R-parity conserving supersymmetry. New R-parity conserving 
SUSY contributions to OuPP decay occur at the level of box diagrams ^^2q- Double beta 
decay then yields information on the mass splitting in the sneutrino-antisneutrino system 
|22j . These constraints leave room for accelerator searches for certain manifestations 
of the second and third generation (B-L)-violating sneutrino mass term, but are most 
probably too tight for first generation (B-L)-violating sneutrino masses to be searched 
for directly. It has been discussed recently jB6j that Oz//5/5 decay by R-parity violating 
SUSY experimentally may not be excluded, although this would require making R-parity 
violating couplings generation dependent. 

We show, in section III that indirect support for the observed evidence for neutrino- 
less double beta decay evidence comes from analysis of other Ge double beta experiments 
(though they are by far less sensitive, they yield independent information on the back- 
ground in the region of the expected signal). 

Table 1: Recent support of the neutrino mass deduced from Oz//3/3 decay P^ElElOIIIl 

by other experiments, and by theoretical work. 



Experiment 


References 


rriy (degenerate i/'s)(eV) 


Ouf3f3 


P El Ellia EH 


0.12 - 0.9 


WMAP 


miEn] 


< 0.23, or 0.33, or 0.50 


CMB 


[IS] 


< 0.7 


CMB+LSS+X-ray gal. Clust. 


[ZH] 


~ 0.2 eV 


SDSS + WMAP 


[HBl 


< 0.57eV 


Z - burst 




0.08 - 1.3 


g-2 




> 0.2 


Tritium 




< 2.2 - 2.8 


u oscillation 


68, 69 


> 0.04 


Theory: 






A4-symmetry 


m 


> 0.2 


identical quark 






and u mixing at GUT scale 




> 0.1 


Alternative cosmological 






'concordance model' 


m 


order of eV 



The discussion in section IV, V, VI, may now just still be of historical interest. Here 
we disprove some criticism of our earlier given results. We show by measurements with 
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a Ra source, performed in 2003 [Tzj, and by various statistical calculations, that the 
criticism by Aalseth et al., (see Mod. Phys. Lett. A17 (2002) 1475-1478), Zdesenko et 
al., (see Phys. Lett. B 546 (2002) 206-215), lanni (in NIM 2004), Feruglio et al., (see 
Nucl. Phys. B 637 (2002) 345) of our earlier results [H 121 E] just was wrong. 

In section VII we give a short discussion, stressing that the evidence for neutrinoless 
double beta decay has been supported by various recent experimental results from other 
fields of research (see Table Q). It is consistent [20] with recent results from cosmic 
microwave background experiments [7^ [TIJ [TH] . The precision of WMAP even allows to 
rule out some old-fashioned nuclear double beta decay matrix elements (see [^). 

It has been shown to be consistent with the neutrino masses required for the Z-burst 
scenarios of high-energy cosmic rays [Z21 IM]- It is consistent with a (g-2) deviating 
from the standard model expectation [^3]. It is consistent also with the limit from the 
tritium decay experiments 03] but the allowed confidence range still extends down to a 
range which cannot be covered by future tritium experiments, if at all jHHI- It is further 
supported by recent theoretical work [7U \ ITU IH^. 

Cosmological experiments like WMAP are now on the level that they can seriously 
contribute to terrestrial research. The fact that WMAP and less strictly also the tritium 
experiments cut away the upper part of the allowed range for the degenerate neutrino 
mass ("^^com ~ 

(0.14 — 3.6) eV) could indicate that the neutrino mass eigenvalues have the 
same CP parity [21 J. 




0.01 



> 0.1 
S 71, 



10 



f lot GENIUS 
I It 



XMASS 

lot 



EXO 

lot 



7kgy W 
0.4 eV — 

" ELEGANT 
BeijinjA 
(76%) [ 



ELEGANT 
(68%) 

UCI TPC 
(68%) 



1 



CUORE 



ot jio' 
A lit 



MOON 

NEMO 3 

2007 I MAMIAND 

10 kg I CaltecH 1 I 
iNeuehalelr I 
Milano' 

TeO^I TPC|o^o[? 
Kiev ' " ' 



1 



UCI 
TPC 



1 



present 
limits 

potential 



potential 
of future 
projects 



48 48 76 76 82 100 100 116 130 136 136 150 
Ca Ca Ge Ge Se Mo Mo Cd Te Xe Xe Nd 



Figure 1: Present sensitivity, and expectation for the future, of the most promising f3f3 
experiments. Given are limits for (m), except for the HEIDELBERG-MOSCOW experi- 
ment where the measured value is given {3a c.l. range and best value). Framed parts of 
the bars: present status; not framed parts: future expectation for running experiments; 
solid and dashed lines: experiments under construction or proposed, respectively. For 
references see [H El Ei • 

Finally we briefly comment in section VIII about the possible future of the field of 
double beta decay. First results from GENIUS-TF which has come into operation on May 
5, 2003 in Gran Sasso with first in world 10 kg of naked Germanium detectors in liquid 
nitrogen jSOlESlEIl, are discussed in another contribution to this report 
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2 Results Obtained in the Period August 2, 1990 
Until May 20, 2003. 

The status of present double beta experiments is shown in Fig. [T] and is extensively 
discussed in The HEIDELBERG-MOSCOW experiment using the largest source 

strength of 11 kg of enriched ^^Ge (enrichment 86%) in form of five HP Ge-detectors is 
running since August 1990 in the Gran-Sasso underground laboratory [121 0111121 1111113111 • 
We present here in Figs. (20 the results obtained with three more years data, until May 
20, 2003. Fig. 121 shows the full spectrum. Fig. IHlthe range around the Qpp value. They 
correspond to a total measuring time of 71.7 kg y. 




1250 1300 1350 1400 1450 1500 1550 1600 

energy [keV] 



2400 2450 2500 2550 2600 2650 2700 2750 2800 
energy [keV] 



Figure 2: The total sum spectrum measured over the full energy range (low-energy part 
(left), and higher energy part (right)) of all five detectors (in total 10.96kg enriched in 
^^Ge to 86%) - for the period 2 August 1990 to 20 May 2003. 

Fig. 121 shows that the line at Qpjs is now - as the Bi lines at 2010.7, 2016.7, 2021.8, 
2052.9 keV - directly clearly seen, while in our first results they had to be projected out 
from the background by a peak search procedure [H 121 ISIHI • 

Earlier measurements of Qpp by [23 EDI EH yielded 2040.71±0.52 keV, 2038.56 ± 
0.32 keV and 2038.668±2.142keV. The precision measurement of [2H1 yields 2039.006 
(50) keV. 

The data have been analysed with various statistical methods. We always process 
background-plus-signal data since the difference between two Poissonian variables does 
not produce a Poissonian distribution [2H]- This is important, but sometimes overlooked 
(see section 6). Analysis of the spectra by nonlinear least squares method, using 
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Figure 3: The total sum spectrum of all five detectors (in total 10.96 kg enriched in ^^Ge), 
for the period August 1990 to May 2003 (71.7kgy) left, and for the period November 
1995-2003 (56.66 kg y) in the range 2000 - 2060 keV and its fit (see section 3.2). 



the Levenberg-Marquardt algorithm yields the fits, shown in Fig. El In these fits the 
peak positions, widths and intensities of all lines are determined simultaneously, and also 
the absolute level of the background. The shape of the latter was chosen to be slightly 
decreasing with energy, corresponding to the complete simulation of the background per- 
formed in by GEANT4. E.g. in Fig. El right, the fitted background corresponds 
to (55.94±3.92) kgy if extrapolated from the background simulated in ^HI for the mea- 
surement with 49.59 kgy of statistics (see Fig. [T^ . This is almost exactly the statistical 
significance of the present experiment (56.66 kg y) and thus a very nice proof of consis- 
tency. Assuming a constant background in the range 2000 - 2060 keV or keeping also the 
slope of a linearly varying background as a free parameter, yields very similar results. 
Analysis with the Maximum Likelihood Method gives results consistent with the above 
method. 

The signal at Q/j^ in the full spectrum (the fit of Fig. El right, yields 2038.44±0.45keV), 
reaches a 4.2o" confidence level for the period 1990-2003, and of 4. la for the period 1995- 
2003 (for details we refer to 0). A detailed description of the analysis of the full data 
1990-2003 will be given in the next Annual Report. 

3 Measurements With a '^^^Bi Source, 

Comparison With Other Ge-Experiments 

By the peak search procedure developped on basis of the Bayes and Maximum Like- 
lihood Methods, exploiting as important input parameters the experimental knowledge 
on the shape and width of lines in the spectrum, weak lines of ^^^Bi had been identified 
at the energies of 2010.7, 2016.7, 2021.6 and 2052.9 keV already in |I1I21 1211101 • Though 
the lines with our improved statistics and analysis are now clearly seen directly in the 
spectrum (Fig. El), we show for comparison the result of the peak search procedure for 
the spectrum taken 1995-2003, in Fig. HI As usual, shown is the probability that there is 
a line of correct width and of Gaussian shape at a given energy, assuming all the rest of 
the spectrum as flat background (which is a highly conservative assumption). 
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Figure 4: Scan for lines in the full spectrum taken from 1995 - 2003 with detectors Nr. 
1,2,3,4,5, with the MLM method (see text). The Bi lines at 2010.7, 2016.7, 2021.8 and 
2052.9 keV are clearly seen, and in addition a signal at ~ 2039 keV. 

Concerning the intensities of these '^^^Bi lines, one has to note that the 2016 keV line, 
as an EO transition, can be seen only by coincident summing of the two successive lines 
E = 1407.98 keV and E = 609.316 keV. Its observation proves that the ^ss-y impurity 
from which it is originating, is located in the Cu cap of the detectors. 

We performed, in the first half of 2003, a measurement of a ^^^Ra source with a high- 
purity germanium detector |T2j. The aim of this work was to investigate the difference 
in the Bi spectra when changing the position of the source with respect to the detector, 
and to verify the effect of TCS (true coincidence summing) for the weak "^^"^Bi lines seen 
in the HEIDELBERG-MOSCOW experiment . 
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Figure 5: Measured ^^^Ra spectrum in the energy range from 2000 to 2100 keV. The upper 
spectrum corresponds to the close geometry, the bottom spectrum to the far geometry. 
The weak lines from ^^^Bi are nicely visible, together with the effect of the true coincidence 
summing at 2016. 7 keV (from [T^l. 

The activity of the ^^^Ra source was 95.2 kBq. The isotope ^^^Ra appears in the 
^^®U natural decay chain and from its decays also ^^^Bi is produced. The 7-spectrum of 
^^^Bi is clearly visible in the ^^^Ra measured spectrum (see Fig. EI). We also performed 
a simulation of our measurement with the GEANT4 simulation tool and we find good 
agreement between the simulation and the measurement ^2]- The premature estimates 



of the Bi intensities given in Aalseth et ai, hep-ex/0202018 and Feruglio et al., Nucl. 
Phys. B 637 (2002), 345, are incorrect, because this long-known spectroscopic effect of 
true coincident summing j2Zj has not been taken into account, and also no simulation of 
the setup has been performed (for details see El El 1121 El)- 

These Bi lines occur also in other investigations of double beta decay. There are three 
other Ge experiments which have looked for double beta decay of ^^Ge. First there is the 
experiment by Caldwell et al. [33], using natural Germanium detectors (7.8% abundance 
of ^^Ge, compared to 86% in the HEIDELBERG-MOSCOW experiment). This was the 
most sensitive natural Ge experiment. With their background a factor of 9 higher than 
in the HEIDELBERG-MOSCOW experiment and their measuring time of 22.6 kg years, 
they had a statistics of the background by a factor of almost four larger than in the 
HEIDELBERG-MOSCOW experiment. This gives useful information on the composition 
of the background. 

Applying the same method of peak search as used in Fig. HI yields (see also [71 ITT]) 
indications for peaks essentially at the same energies as in Fig. [3] (see Fig. This shows 
that these peaks are not fluctuations. In particular it sees the 2010.78, 2016.7, 2021.6 and 
2052.94 keV ^^^Bi lines, but also the unattributed lines at higher energies. It finds, 
however, n o line at 2039 keV. This is consistent with the expectation from the rate found 
in the HEIDELBERG-MOSCOW experiment. About 29 identified events observed during 
1990-2003 in the latter correspond to 0.7 expected events in the Caldwell experiment, 
because of the use of non-enriched material and the shorter measuring time. Fit of the 
Caldwell spectrum allowing for the ^^"^Bi lines and a 2039 keV line yields 0.4 events for the 
latter (see [3] and Fig. [HJ). 




Figure 6: Result of the peak-search procedure performed for the UCBS/LBL spectrum 
(left: Maximum Likelihood method, right: Bayes method). On the y axis the proba- 
bility of having a line at the corresponding energy in the spectrum is shown (from [T] ) • 

The first experiment using enriched (but not high-purity) Germanium 76 detectors 
was that of Kirpichnikov and coworkers [HHl- These authors show only the energy range 
between 2020 and 2064 keV of their measured spectrum. The peak search procedure 
finds also here indications of lines around 2028 keV and 2052 keV (see Fig. but not 
any indication of a line at 2039 keV. This is consistent with the expectation, because for 
their low statistics of 2.95 kgy they would expect here (according to HEIDELBERG- 
MOSCOW) 1.1 counts. 
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Figure 7: Result of the peak-search procedure performed for the ITEP/YePI spectrum 
[SS! (from 0111]). 
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Figure 8: Result of the peak-search procedure performed for the IGEX spectrum [[57] . 
Left: Maximum Likelihood method, right: Bayes method. On the y axis the probability 
of having a line at the corresponding energy in the specrtum is shown (from [71 ITT]). 
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Figure 9: Analysis of the spectrum measured by D. Caldwell et al. [21], with the 
Maximum Likelihood Method, in the energy range 2000-2060 keV assuming lines at 2010.7, 
2016.7, 2021.6, 2052.9, 2039.0 keV. No indication for a signal at 2039 keV is observed in 
this case (see (Tj). 
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Another experiment (IGEX) used between 6 and 8.8 kg of enriched ''^Ge, but collected 
since beginning of the experiment in the early nineties till shutdown in 1999 only 8.8 
kg years of statistics jSTj. The authors of unfortunately show only the range 2020 
to 2060 keV of their measured spectrum in detail. Fig. |H1 shows the result of our peak 
scanning of this range. Clear indications are seen for the lines at 2021 and 2052 keV, but 
also of the unidentified structure around 2030 keV. Because of the conservative assumption 
on the background treatment in the scanning procedure (see above) there is no chance 
to see a signal at 2039 keV because of the 'hole' in the background of that spectrum (see 
Fig. 1 in jHZ!)- With some good will one might see, however, an indication of ~3 events 
here, consistent with the expectation of the HEIDELBERG-MOSCOW experiment 
of ~ 2.6 counts. 

4 Statistical Features: 

Sensitivity of Peak Search, Analysis Window 

For historical reasons, at this point it may be useful to demonstrate the potential of 
the peak search procedure used in [U El El- Fig. ITUl shows a spectrum with Poisson- 
generated background of 4 events per channel and a Gaussian line with width (standard 
deviation) of 4 channels centered at channel 50, with intensity of 10 (left) and 100 (right) 
events, respectively. Fig. ^1 shows the result of the analysis of spectra of different line 
intensity with the Bayes method (here Bayes 1-4 correspond to different choice of the prior 
distribution: (1) fi{ri) = 1 (flat), (2) fi{ri) = 1/1], (3) fi{ri) = l/y/fj, (4) Jeffrey's prior) and 
the Maximum Likelihood Method. For each prior 1000 spectra have been generated with 
equal background and equal line intensity using random number generators available at 
CERN [211. The average values of the best values agree (see Fig. IT^ very well with the 
known intensities also for very low count rates (as in Fig. ^1 left). 




channel channel 

Figure 10: Example of a random-generated spectrum with a Poisson distributed background 
with 4.0 events per channel and a Gaussian line centered in channel 50 (line-width corresponds 
to a standard-deviation of cr = 4.0 channels). The left picture shows a spectrum with a line- 
intensity of 10 events, the right spectrum a spectrum with a line-intensity of 100 events. The 
background is shown dark, the events of the line bright (see jllj). 

In Fig. ^Iwe show two simulations of a Gaussian line of 15 events, centered at channel 
50, again with width (standard deviation) of 4 channels, on a Poisson-distributed back- 
ground with 0.5 events/channel. The figure gives an indication of the possible degree of 
deviation of the energy of the peak maximum from the transition energy, on the level of 
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Figure 11: Result of an analysis as function of the evaluation width. The used spectrum 
consists of a Poisson distributed background with 4 events per channel, and a line of 10 
events (see Fig. ^1 left part). The dark area corresponds to a 68.3% confidence area with 
the dark line being the best value. Below an evaluation width of 35 channels the result 
becomes unreliable, above 35 channels the result is stable (see also [71ITT]). 



statistics collected in experiments like the HEIDELBERG-MOSCOW experiment (here 
one channel corresponds to 0.36 keV). This should be kept in mind. 
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Figure 12: Results of analysis of random- number generated spectra, using Bayes and Maximum 
Likelihood method (the first one with different prior distributions). For each number of events 
in the simulated line, shown on the x-axis, 1000 random generated spectra were evaluated with 
the five given methods. The analysis on the left side was performed with an Poisson distributed 
background of 0.5 events per channel, the background for the spectra on the right side was 
4.0 events per channel. Each vertical line shows the mean value of the calculated best values 
(thick points) with the la error area. The mean values are in good agreement with the expected 
values (horizontal black dashed lines) (see j3 ^| ) . 

The influence of the choice of the energy range of the analysis around Q^js has been 
thoroughly discussed in [21 E] ■ Since erroneous ideas about this point are still around, 
let us remind of the analysis given in [21 El U\ which showed that a reliable result is 
obtained for a range of analysis of not smaller than 35 channels (i.e. ±18 channels) - one 
channel corresponding to 0.36 keV in the HEIDELBERG-MOSCOW experiment (see Fig. 
[TT|) . This is an important result, since it is, in case of a weak signal, of course important to 
keep the range of analysis as s m a 1 1 as possible, to avoid to include lines in the vicinity 
of the weak signal into the background (see, e.g. Fig. 9 in jHI])- This unavoidably occurs 
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Figure 13: Two spectra with a Poisson-distributed background and a Gaussian line with 15 
events centered in channel 50 (with a width (standard-deviation) of 4.0 channels) created with 
different random numbers. Shown is the result of the peak-scanning of the spectra. In the left 
picture the maximum of the probability corresponds well with the expected value (black line) 
whereas in the right picture a larger deviation is found. When a channel corresponds to 0.36 keV 
the deviation in the right picture is ~ 1.44 keV (see [T} 111)). 



when e.g. proceeding as suggested in F. Feruglio et al., hep-ph/0201291 and Nucl. Phys. 
B 637 (2002) 345-377, Aalseth et. al., hep-ex/02020T8'and Mod. Phys. Lett. A 17 (2002) 
1475, Yu.G. Zdesenko et. al., Phys. Lett. B 546 (2002) 206, A. lanni, in Press NIM 
2004. The arguments given in those papers are therefore incorrect. Also Kirpichnikov, 
who states that his analysis finds a 2039 keV signal in the HEIDELBERG-MOSCOW 
spectrum on a 4 sigma confidence level (as we also see it) makes this mistake, when 
analysing the pulse shape spectrum. 

The above discussion is now in this context only of historical interest, since with the 
better statistics we have now, we can analyze simultaneously a large energy range (as 
shown in Fig. EJ. 



5 Simulation with GEANT4 



Finally the background around (5/3/3 will be discussed from the side of simulation. A 
very careful new simulation of the different components of radioactive background in the 
HEIDELBERG-MOSCOW experiment has been performed by a new Monte Carlo pro- 
gram based on GEANT4 ^3] ■ This simulation uses a new event generator for simulation 
of radioactive decays basing on ENSDF-data and describes the decay of arbitrary radioac- 
tive isotopes including alpha, beta and gamma emission as well as conversion electrons 
and X-ray emission. Also included in the simulation is the influence of neutrons in the 
energy range from thermal to high energies up to 100 MeV on the measured spectrum. 
Elastic and inelastic reactions, and capture have been taken into account, and the corre- 
sponding production of radioactive isotopes in the setup. The neutron fluxes and energy 
distributions were taken from published measurements performed in the Gran Sasso. Also 
simulated was the cosmic muon flux measured in the Gran Sasso, on the measured spec- 
trum. To give a feeling for the quality of the simulation. Fig. IT^ shows the simulated and 
the measured spectra for a '^'^^Th source spectrum for as example one of our five detectors. 
The agreement is excellent. 
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Figure 14: Comparison of the measured data (black line, November 1995 to April 2002) 
and simulated spectrum (red line) for the detectors Nrs. 1,2,3 and 5 for a ^^^Th source 
spectrum. The agreement of simulation and measurement is excellent (from [T5]). 

The simulation of the background of the experiment reproduces a 1 1 lines observed 
in the sum spectrum of the five detectors, in the energy range between threshold (around 
100 keV) and 2020 keV Hg. 




E [keV] 

Figure 15: Simulated background of the HEIDELBERG-MOSCOW experiment in the 
energy range from 2000 to 2100 keV with all known background components, for the 
period 20 November 1995 to 16 April 2002 (from [13]). 

Fig. fTHl shows the simulated background in the range 2000-2100 keV with all 
known background components. 

The background around is according to the simulations flat, the only expected 
lines come from ^^^Bi (from the ^^^U natural decay chain) at 2010.89, 2016.7, 2021.6, 
2052.94, 2085.1 and 2089.7keV. Lines from cosmogenically produced ^*^Co (at 2034.76 keV 
and 2041.16 keV), half-life 77.3 days, are not expected since the first 200 days of measure- 
ment of each detector are not used in the data analysis. Also the potential contribution 
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from decays of ''^Ge, ^^Ga, or ^^^Ac, should not lead to signals visible in our measured 
spectrum near the signal at For details we refer to jlSI- 

6 Proofs and disproofs 

Our earlier result published in P3 12113) which now is confirmed on a Aa level, had been 
questioned in some papers [Aalseth et al, hep-ex/0202018, and in Mod. Phys. Lett. A 
17 1475-1478; Feruglio et al, Nucl. Phys. B 637 (2002) 345; Zdesenko et al., Phys. 
Lett. B 546 (2002) 206], and Kirpichnikov, talk at Meeting of Physical Section of Russian 
Academy of Sciences, Moscow, December 2, 2002, (and priv. communication, Dec. 3, 
2002) and A. lanni, Nucl. Instruments A (2004) (available online 28 September 2003). 
We think that we have shown in a convincing way during 2002 and 2003 that these claims 
against our results were incorrect in various ways, and have published our arguments in 
[T^ ITU [71 . In particular the estimates of the intensities of the lines in the first 

two papers do not take into account the effect of true coincidence summing, which can 
lead to drastic underestimation of the intensities. A correct estimate would also require a 
Monte Carlo simulation of our setup, which has not been performed in the above papers. 

The paper by Zdesenko et al. starts from an arbitrary assumption, namely that 
there are lines in the spectrum at best only at 2010 and 2053 keV. This contradicts to 
the experimental result, according to which there are further lines in the spectrum (see 
Fig. El in this report). For example they could have easily deduced from the intensity of 
the 2204 keV Bi line in the measured spectrum (Fig. 12}) that lines at 2053 keV etc. are 
expected [32^ . In this way and also by some subtraction procedure, ignoring that the result 
of subtracting a Poisson-distributed spectrum from a Poisson-distributed spectrum does 
not give a Poisson distributed spectrum (see, e.g. |nS|) they come to wrong conclusions. 

Kirpichnikov states [36j that from his analysis he clearly sees the 2039 keV line in the 
full (not pulse-shape discriminated) spectrum on a Aa level. He claims that he does not 
see the signal in the pulse shape spectrum. The simple reason to see less intensity is that 
in this case he averages for determination of the background over the full energy range 
without allowing for any lines. 

All of these papers, when discussing our earlier choice of the width of the search window 
(in the analysis of the data taken until May 2000), ignore the results of the statistical 
simulations - we present here, and have published in i2|i2l lH El ITUl ITT|. 

The strange effects found recently by the Kurchatov people [221 in their rough analysis 
of part of the data, have been traced back to including corrupt data into the analysis. 
The artefacts seen in their Figs. 4,5,7,8 do not exist in our data, which lead to the results 
shown in Figs. I2I3I ffor details see [HI I3E]). 

7 Discussion of results 

We emphasize that we find in all analyses of our spectra a line at the value of Q/3/3. The 
results confirm our earlier result with higher statistics. For details we refer to the next 
Annual Report and to 6]. 

The result obtained is consistent with all other double beta experiments - which 
reach in general by far less sensitivity. The most sensitive experiments following the 
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HEIDELBERG-MOSCOW experiment are the geochemical ^^^Te experiment with ^ 
2(7.7) X 10^4 y (68% c.L), EH the ^^^Xe experiment by the DAMA group with T°y2 > 
1.2 X 10^^ y (90% C.I.), a second enriched ^^Ge experiment with T°^2 > 1-2 x lO^^ y ^ and 
a "°*Ge experiment with Tj/g > 1 x 10^^ y 34 . Other experiments are aheady about a 
factor of 100 less sensitive concerning the 0i>/3f3 half-hfe: the Gotthard TPC experiment 
with i^^Xe yields T°j2 > 4-4 x lO^^ y (90% c.l.) and the Milano Mibeta cryodetector 
experiment T°j2 > 1-44 x 10^^ y (90% c.l.). 

Another experiment '^ET] with enriched ^^Ge, which has stopped operation in 1999 after 
reaching a significance of 8.8 kgy, yields (if one believes their method of 'visual inspection' 
in their data analysis), in an analysis correcting for on arithmetic error which has been 
made in ^ (for discussion see [SB]) a hmit of about T°j2 > 5 x 10^^ y (90% c.l.). The 
^28Te geochemical experiment yields (m^) < 1.1 eV (68 % c.l.) 071, the DAMA ^^^Xe 
experiment (mj,) < (1.1—2.9) eV and the ^^°Te cryogenic experiment yields (mjy) < 1.8 eV. 

Concluding we obtain, with > Aa probability, evidence for a neutrinoless double beta 
decay signal. Following this interpretation, at this confidence level, lepton number is 
not conserved. Further the neutrino is a Majorana particle. If the Oz//3/3 amplitude is 
dominated by exchange of a massive neutrino the effective mass (m) is deduced from the 
full spectrum (using the matrix elements of [23]) to be (m) = (0.1 - 0.9) eV {3a confidence 
range), allowing already for a ± 50% uncertainty of the matrix element. The best value 
is 0.4 eV. 

Assuming other mechanisms to dominate the Ou/SlS decay amplitude, the result allows 
to set stringent limits on parameters of SUSY models, leptoquarks, compositeness, masses 
of heavy neutrinos, the right-handed W boson and possible violation of Lorentz invariance 
and equivalence principle in the neutrino sector. For a discussion and for references we 
refer to [H] |S3 dHl ISS HE] • 

With the value deduced for the effective neutrino mass, the HEIDELBERG-MOSCOW 
experiment excludes several of the neutrino mass scenarios allowed from present neutrino 
oscillation experiments (see Fig. [TB|) - allowing only for a degenerate mass scenario fSl 
l2nilS]- Fig- IT^ shows also the limits obtained from WMAP, which at the present level of 
sensitivity is not able to rule out any neutrino mass scheme. 

The evidence for neutrinoless double beta decay has been supported by various recent 
experimental and theoretical results (see Table 1). Assuming the degenerate scenarios 
to be realized in nature we fix - according to the formulae derived in pHI - the common 
mass eigenvalue of the degenerate neutrinos to m = (0.1 - 3.6) eV. Part of the upper 
range is excluded by tritium experiments, which give a limit of m < (2.2 - 2.8) eV (95% 
c.l.) [13]. The full range can only partly (down to ~ 0.5 eV) be checked by future 
tritium decay experiments, but might be checked by some future (3/3 experiments (see 
next section). Recent theoretical work [HO] even doubts, that tritium experiments are in 
principle capable to check a Oz//?/3 result. The deduced best value for the mass is consistent 
with expectations from experimental /i — 67 branching limits in models assuming the 
generating mechanism for the neutrino mass to be also responsible for the recent indication 
for as anomalous magnetic moment of the muon [65jj . It lies in a range of interest also for 
Z-burst models recently discussed as explanation for super-high energy cosmic ray events 
beyond the GKZ-cutoff [Mj 17^ and requiring neutrino masses in the range (0.08 - 1.3) 
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Figure 16: The impact of the evidence obtained at {4.2a c.L) for neutrinoless double beta 
decay (best value of the effective neutrino mass (m) = 0.4 eV, 3a confidence range (0.1 
- 0.9) eV - allowing already for an uncertainty of the nuclear matrix element of a factor 
of ± 50%) on possible neutrino mass schemes. The bars denote allowed ranges of (m) 
in different neutrino mass scenarios, still allowed by neutrino oscillation experiments (see 
[13 I2ni)- AH models except the degenerate one are excluded by the new Oz//3/3 decay 
result. Also shown is the exclusion line from WMAP, plotted for J2'^u < 1.0 j76j . 
WMAP does not rule out any of the neutrino mass schemes. Further shown are the 
expected sensitivities for the future potential double beta experiments CUORE, MOON, 
EXO and the 1 ton and 10 ton project of GENIUS [HI UHl EOl iH (from [20]). 

eV. A recent model with underlying A4 symmetry for the neutrino mixing matrix also 
leads to degenerate Majorana neutrino masses > 0.2 eV, consistent with the present result 
from OuPP decay [nZHZOj. The result is further consistent with the theoretical paper of 
j71j . Starting with the hypothesis that quark and lepton mixing are identical at or near 
the GUT scale, Mohapatra et al. [21 show that the large solar and atmospheric neutrino 
mixing angles can be understood purely as result of renormalization group evolution, 
if neutrino masses are quasi-degenerate (with same CP parity). The common Majorana 
neutrino mass then must be, in this model, larger than 0.1 eV. An completely independent 
theoretical proof, that neutrinos should have Majorana nature, has been given recently 



For WMAP a limit for the sum of the neutrino masses of = < 0.69 eV at 95% 
c.L, was given by the analysis of ref. [71|. More realistically this limit on the total mass 
should be jTB] nis = Y^rrii < 1.0 eVat 95%c.l. The latter analysis also shows, that four 
generations of neutrinos are still allowed and in the case of four generations the limit on 
the total mass is increased to 1.38 eV. If there is a fourth neutrino with very small mass, 
then the limit on the total mass of the three neutrinos is even further weakened and there 
is essentially no constraint on the neutrino masses. In our Fig. ^Iwe show the contour 



by [H21- 
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line for WMAP assuming < 1-0 eV. 

A recent analysis of the Sloan Digital Sky Survey, together with WMAP yields jSS] 

m^ = ^mi<1.7eV at 2 a. (8) 

Comparison of the WMAP results with the effective mass from double beta decay rules out 
completely (see jZH]) a 15 years old old-fashioned nuclear matrix element of double beta 
decay, used in a recent analysis of WMAP [77j . In that calculation of the nuclear matrix 
element there was not included a realistic nucleon-nucleon interaction, which has been 
included by all other calculations of the nuclear matrix elements over the last 15 years. 

The range of (m) fixed in this work is, in the range to be explored by the satellite 
experiments MAP and PLANCK [TH 17^ 17^. The limitations of the information from 
WMAP are seen in Fig. Ilfif thus results of PLANCK gerly awaited. 

The neutrino mass deduced leads to 0.002> Q^h"^ < 0.1 and thus may allow neutrinos 
to still play an important role as hot dark matter in the Universe |47j . 

8 Future of P/S experiments 

With the HEIDELBERG-MOSCOW experiment, the era of the small smart experiments 
is over. New approaches and considerably enlarged experiments (as discussed, e.g. in 
daEnmniliainniillllliZl) wm be required in future to fix the Oz//3/3 half life of 
^^Ge with higher accuracy. This will, however, because of the uncertainties in the nuclear 
matrix elements, which probably hardly can be reduced to less than 50%, only marginally 
reduce the precision of the deduced neutrino mass. 

Since it was realized in the HEIDELBERG-MOSCOW experiment, that the remaining 
small background is coming from the material close to the detector (holder, copper cap, 
...), elimination of any material close to the detector will be decisive. Experiments which 
do not take this into account, will allow at best only rather limited steps in sensitivity. 
Furthermore there is the problem in cryodetectors that they cannot differentiate between 
a P and a 7 signal, as this is possible in Ge experiments. 

Another crucial point is the energy resolution, which can be optimized only in exper- 
iments using Germanium detectors, or, to some less extent, with bolometers. It will be 
difficult to probe evidence for this rare decay mode in experiments, which have to work - 
as result of their limited resolution - with energy windows around Q/3/3 of several hundreds 
of keV. 

Another important point is the efficiency of a detector for detection of a signal. 
For example, with 14% efficiency a potential future 100 kg ^"^Se experiment would be, 
because of its low efficiency, equivalent only to a 10 kg experiment (not talking about the 
energy resolution). 

In the first proposal for a third generation double beta experiment, our GENIUS 
proposal [Sni 113 HOI 1121 inni El , the idea is to use 'naked' Germanium detectors in a huge 
tank of liquid nitrogen. It seems to be at present the only proposal, which can fulfill 
both requirements mentioned above - to increase the detector mass and simultaneously 
reduce the background drastically. At the present status of results of the HEIDELBERG- 
MOSCOW experiment, however - with a confidence level of ~ 4cr, it is questionable, 
whether GENIUS would be needed for decay. Probably it would be preferable to 
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perform an experiment with another isotope but fulfilling all requirements mentioned 
above. The GENIUS-Test-Facility, originally planned to prove the feasibility of some key 
constructional parameters of GENIUS, and put into operation on May 5, 2003 in GRAN 
SASSO, could however, play an important role in testing the evidence seen for cold 
dark matter by DAMA (see |^ l52j. and another Report to this volume). Only a GENIUS 
with some ten tons of enriched ^^Ge might possibly be of interest, to investigate whether 
another exotic mechanism such as exchange of SUSY particles, (see, e.g. fH]) niight 
contribute to the Ou^P decay amplitude. This may be, however, a very far dream. 

9 Summary 

The HEIDELBERG-MOSCOW experiment has been continued regularly in 2003. Un- 
fortunately, it had to stop operation according to non-prolongation of our contract with 
Kurchatov institute, at 30 November 2003. Since then still various calibration measure- 
ments with radioactive sources are going on. 

The first analysis of the full data taken with the HEIDELBERG-MOSCOW experi- 
ment in the period 2 August 1990 until 20 May 2003 is presented. The improved statistics 
and data analysis leads to a ~ 4cr evidence for a signal at the Q-value for neutrinoless 
double beta decay. This confirms our earlier claim [TJ |21 13 E]- Additional support for 
this evidence has been presented by showing consistency of the result - for the signal, 
and for the background - with other double beta decay experiments using non-enriched 
or enriched Germanium detectors (see also jZlE])- In particular it has been shown that 
the lines seen in the vicinity of the signal are seen also in the other experiments. This is 
important for the correct treatment of the background. Furthermore, the sensitivity of the 
peak identification procedures has been demonstrated by extensive statistical simulations. 
It has been further shown by new extensive simulations of the expected background by 
GEANT4, that the background around Qjsp should be flat, and that no known gamma 
line is expected at the energy of (see [12] )■ The 2039 keV signal is seen o n 1 y in the 
HEIDELBERG-MOSCOW experiment, which has a by far /arg'er statistics than all other 
double beta experiments. 

The importance of first evidence for violation of lepton number and of the Majorana 
nature of neutrinos is obvious. It requires beyond Standard Model Physics on one side, 
and may open a new era in space-time structure [021 ■ It has been discussed that the 
Majorana nature of the neutrino tells us that spacetime does realize a construct that is 
central to construction of supersymmetric theories. 

One of the consequences of the result of the HEIDELBERG-MOSCOW experiment on 
the present confidence level, may be, that to obtain deeper information on the process of 
neutrinoless double beta decay, new experimental approaches, different from all, what is at 
present persued, may be required. The unique importance of double beta decay to inves- 
tigate the neutrino mass, is stressed by the recent observation, that tritium experiments 
might suffer from principle problems to see a neutrino mass at all [HO] • 

With the successful start of operation of GENIUS-TF with the first four naked Ge 
detectors in hquid nitrogen on May 5, 2003 in GRAN SASSO, which is described in |i4^ l3U] 
(see our second contribution to this Report) a historical step has been achieved of a novel 
technique and into a new domain of background reduction in underground physics in the 
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search for rare events. In the hght of the above comments, natural task of GENIUS-TF 
will be to look for cold dark matter by the modulation signal. 
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